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Abstract: The properties of the Gabor and Morlet transforms are examined with respect to 
the Fourier analysis of discretely sampled data. Forward and inverse transform pairs based 
on a fixed window with uniform sampling of the frequency axis can satisfy numerically the 
energy and reconstruction theorems; however, transform pairs based on a variable window 
or nonuniform frequency sampling in general do not. Instead of selecting the shape of the 
window as some function of the central frequency, we propose constructing a single window 
with unit energy from an arbitrary set of windows which is applied over the entire frequency 
axis. By virtue of using a fixed window with uniform frequency sampling, such a transform 
satisfies the energy and reconstruction theorems. The shape of the window can be tailored 
to meet the requirements of the investigator in terms of time/frequency resolution. The 
algorithm extends naturally to the case of nonuniform signal sampling without modification 
beyond identification of the Nyquist interval. 

Keywords: Fourier transform, Gabor transform, Morlet transform, multiresolution analysis 



1. Introduction 

The primary criticism leveled at the use of the continuous wavelet transform for the spectral analysis 
of discretely sampled data is that it fails to give quantitatively meaningful results. The power spectral 
density produced from the convolution of a wavelet basis and a discrete signal gives a qualitative picture 
of the temporal variation of its frequency content; however, when one attempts the reconstruction of the 
signal, the residual is not on the order of the precision of one's computational device. Likewise, the 
integrated margins of the power spectral density do not precisely equal the energy of the original signal. 
In other words, the discrete implementation of the continuous wavelet transform and its inverse do not 
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satisfy the energy and reconstruction theorems of spectral analysis. The goal of this investigation is to 
devise a multiresolution analysis which does satisfy those theorems. 

Our insistence upon the satisfaction of the energy and reconstruction theorems is because they are first 
principles requirements related to the conservation of physical energy. When suitably defined, energy 
is not to be created nor destroyed; neither should it get lost in the shuffle. There is a deep relation 
between energy content and information content, as quantum mechanics teaches us, thus a loose grip on 
one implies a loose grip on the other. On a more practical level, the satisfaction of those theorems is 
among the requirements for a maximum entropy spectral analysis of data which includes the effects of 
measurement uncertainty, whose consideration is beyond the scope of this article. 

The review by Torrence and Compo [1] remains a popular resource for practitioners of wavelet 
analysis. It relies on the method by Farge [2] for the reconstruction of the data signal. Grossmann 
and Morlet [3] are credited with establishing the reconstruction theorem in the continuum. Meyer [4] 
and Mallat [5] developed the theory of multiresolution analysis, and Daubechies [6] constructed the first 
orthonormal basis with compact support, leading to the implementation of the dyadic wavelet transform 
in terms of finite impulse response digital filter banks [7]. Some interesting applications of the continuous 
wavelet transform for the spectral analysis of data can be found in references [8-19]. 

This paper is organized as follows. First, we will quickly review the theory of the continuous Fourier 
transform and its discrete implementation. Next, we will look at the Gabor transform and its relation to 
the wavelet transform using the Morlet basis. We will then propose an algorithm for a layered window 
transform whose spectral density is similar to that of the Morlet transform yet which satisfies the energy 
and reconstruction theorems. Following that, we will demonstrate that the algorithm works unaltered 
for data with irregular sampling once the corresponding Nyquist interval is identified. Finally, we will 
look at how the selection of the window affects the time/frequency resolution of the transform. We will 
conclude with a brief summary and suggestions for applications. 

Some notations used throughout this paper are explained here. Scalars are written as s, while vectors 
and matrices are written as v and M respectively and may be defined in terms of their components, e.g. 
v = v s = v(s). The transpose of a vector or matrix is indicated by the superscript v T , and the conjugate 
transpose by M^, using the standard rules for matrix multiplication. Inner products may be written in 
bra-ket notation as (a\b) s , and matrix entries as the ket-bra \a)(b\. Expectation values are notated as 
(v(s)) . Sets will be indicated by their boundaries [a, b], and whether continuous or discrete must be 
derived from context. The operation of rounding up, i.e. taking the next greatest integer, is denoted by 
\a] , and when the notation a — > b is encountered, it is understood to mean "a is replaced by b". The 
programs used for this analysis are available from http://www.alphawaveresearch.com online. 

2. Continuous and Discrete Fourier Transforms 

Let us begin by considering the Fourier transform in the continuum over axes of time and frequency. 
Suppose we have some signal y(t), possibly complex valued, of finite energy E y = (y\y) t = 
JZ^o y*(t)y(t)dt < oo. If the signal carries units of u y , then the signal energy has units of ue = u y u t . 
The units of signal energy are proportional to those of physical energy in joules by a factor of the load 
impedance Ej = E y jZ L \ for example, if the signal has units of volts u y = V and time is measured 
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in seconds u t = s so that ue = V 2 s, then accounting for the load impedance in ohms uz = ^ gives 
physical units of V 2 s/Vt = (ml 2 /t 2 q) 2 t/(ml 2 /tq 2 ) in terms of mass, length, time, and charge, which 
is equivalent to joules. The units of frequency are reciprocal to those of time Uf = u^ 1 , and the term 
"power" is understood to mean "energy density" and must be qualified by the domain for its distribution. 
The squared modulus |y(t)| 2 = P(t) can thus be identified as the temporal power of the signal. 
Under suitable conditions not elaborated here, the Fourier integral and its inverse, 

/oo /»oo 
exp(-i2nft)y(t)dt , y(t) = / exp(i27rft)y(f)df , (1) 

define an integral transform pair satisfying the Plancherel energy theorem \y(t)\ 2 dt = \y{f)\ 2 df 
and the Fourier inversion theorem y(t) = y(t). The representation over the frequency axis y(f) carries 
units of Uy = u y u t so that the squared modulus \y(f)\ 2 = P(f) gives the spectral power of the signal, 
and the spectral energy is E y = (y\y)f = y*(f)y{f)df. When the signal is of finite duration 
t G [0, T], the Fourier transform's response to a component sinusoid y(t) = exp(i2irf't) is 

y(f) = f eM-*Mf ~ f'Wt = exp(-m/ A T) Sm(7r/AT) , (2) 
Jo ^/a 

expressed in terms of the frequency offset f A = f — f, thus the spectral power of the uniform window 
function $t(^) = T -1 / 2 normalized to unit energy describes the leakage in the frequency spectrum 
induced by the finite duration of the signal, 

2 

df A , (3) 



L T (fA)df A = \^ T (fA)\ 2 df A =T- 1 



sin(7r/ A T) 



with domain f A G [—00,00]. The continuum leakage function is not periodic in f A , as its magnitude 
decays according to f A 2 . 

With regard to practical data analysis, let us suppose now that the signal is given in terms of discrete 
samples in time with uniform duration and spacing. (Those requirements will be relaxed in Section 5.) 
The time axis can then be described in terms of integers t G [1, T) with unit u t = A t given by the sample 
rate, allowing the signal to be written as a vector y = y t = y(t). The signal energy can be expressed in 
terms of matrix multiplication as E y = y^D 4 y, where the temporal metric D t = l T Aj is proportional to 
the identity matrix of order T. The effect of frequency aliasing induced by the uniform sampling is often 
misunderstood. If we evaluate the Fourier transform of the discrete window with unit energy, 



M/a) = T ' 1/2 E ex P H2vr/ A t) A t = T' 1 ' 2 ex P [-m/ A (T + 1)] ^ f f^ A t 
^ sm(7r/ A J 



t=i 

we find that the leakage function becomes periodic in f A , 

"sin(7r/ A T) 



(4) 



M/a)g?/a ^ T- 1 



_ sm(nf A ) _ 



A 2 df A , (5) 



with period 1 in units of Uf = A t 1 . The principle branch is usually chosen as / G [—1/2, l/2]A t 1 , but 
one should always respect the periodic nature of the spectrum. The normalization of the spectral power 
likewise is defined over one period of the frequency axis, J2t=i \v(t) I — j_i/ 2 Ivif) I In Figure 1 
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Figure 1. Comparison of the continuous time leakage function in (a) to its discrete time 
counterpart in (b) for a uniform window of duration T = 5 over a continuous frequency axis. 





we compare the continuum and discrete leakage functions for a window with duration T = 5, where we 
see the expected oscillatory lobes as well as the periodicity induced by the discrete sampling in time. 

Now let us consider the discretization of the frequency axis into uniform bins over its principle branch 
/ E [— fcfc), where f c = (2Ai) _1 is the Nyquist critical frequency. To fully specify the frequency 
metric D/, one must state its order N, equal to the number of positive frequencies < f c , as well as its 
parity P E {0, 1}, indicating whether an even or odd number of bins N' = 2N + P is used to span 
the domain. When P = 1, the edgemost bins corresponding to frequencies ±/ c have a bin width equal 
to one half that of the others, so that the limits of integration are respected; because of aliasing, the 
integrand will have the same value at those locations, so that only one full bin's contribution is counted. 
This convention differs from that usually given for the discrete Fourier transform [20], where the factor 
of 1/2 is absorbed by the coefficients rather than the metric. The order N specifies the spacing of the 
frequencies as A/ = (2NA t y 1 , and the metric can be written as D f /A f = ljv'-(| | + |iV / )(A^ / |)P/2 
so that Tr D/ = A^ 1 . The central frequencies of the bins can be expressed as f n =[n — (N' + l)/2] A/ 
for integers n E [1, N'}. Figure 2 compares the even and odd discretizations of order N = A. While the 
odd parity discretization is perhaps more familiar, for many of our purposes the even discretization will 
prove more convenient. 

We are finally ready to define the two-sided discrete Fourier transform (of order N and parity P) of 
a discretely sampled signal (of duration T). If we collect our Fourier basis functions into the form of 
a matrix, = 0(t, n) = ex.p(i2nf n t), then we can easily write the discrete transform pair in terms 
of matrix multiplication, y = Q^D t y and y = 0D^y, as well as the spectral energy E$ = y^D^y. In 
component notation, one has 



N' 



y{n) 



^exp(-z27r/ n t)y(£)A t , y(t) = ^ exp(i2nf n t)y(n)A f 



(6) 



71=1 



where Aj is understood to account for the edge bins if P = 1. If the signal y{t) is real, so that 
y(—f) = y(f)*, then one may define the one-sided transform, retaining only the nonnegative portion of 
the frequency axis with N' — >• N + P bins and /„ — ► [n — (P + l)/2]Aj, by renormalizing the basis 
functions to twice the energy via — > y/2Q and letting y(t) — > Re y(t). Hereafter, we will assume y(t) 
is real so that we can focus our attention on the one-sided spectrum. 
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Figure 2. Comparison of the even (a) and odd (b) frequency discretizations of order N = 4 
with the central frequencies indicated by x . 

N = 4, P = N = 4, P = 1 
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In the fully discrete setting, the satisfaction of the energy and reconstruction theorems is dependent 
upon there being a sufficient number of degrees of freedom (dof) in the basis functions to fully represent 
the information content of the signal. For real y(t), the number of dof is equal to the signal duration 
T. For the one-sided transform, whether of even or odd parity, the number of dof is equal to 2N, as 
each frequency's coefficient has both an amplitude and a phase, y n = A n exp(iuj n ), except for the cases 
f n = or f c which have only an amplitude. The critical frequency f c is identified as the lowest positive 
frequency whose basis function is entirely real over the discrete time axis, an observation which will 
prove useful when we consider the case of irregularly sampled data. The minimal order beyond which 
the energy and reconstruction theorems are satisfied can be evaluated as A^ min = \T/2] for a real signal 
(and N min = T for a complex one). That condition is realized (for even T) when 0^0 A t = \n'T, 
indicating that the basis functions form an orthogonal set for N = T/2. However, orthogonality over 
the discrete time axis is not a requirement, as y^D t y = y^D^y and y = y exactly (meaning to the 
precision of one's computational device) for any N > N miQ and either parity P. 

Lastly, let us look at what happens when one tries to express the fully discrete, one-sided Fourier 
transform over an axis of period (or scale) s rather than frequency /. In the continuum, the relation 
between the axes / = s -1 yields the Jacobian \df/ds\ = s~ 2 , thus the spectral power over scale P(s) = 
\y{s)\ 2 is equal to the spectral power over frequency P(f) multiplied by the Jacobian, P(s) = f 2 P(f), 
otherwise expressed as y(s) = fy(f). In the discrete setting, however, one must be explicit with the 
mapping of the bin boundaries so that J2 n \y(f n )\ 2 Af = J2 n \y(s n )\ 2 A s . Let A/ = b — a be the width of 
one frequency bin centered on f n — (a + b)/2, which gets mapped to the scale bin A s = a -1 — b^ 1 such 
that Aj/A s = ab and P(s n ) = abP(f n ). One may very well ask, in what sense is s n = f~ l the center of 
the scale bin? The value f n is both the mean and the median of the frequency bin with uniform measure 
AT 1 , but the mean period over that bin is {f~ l ) f = log(6/a) AT 1 . The value s n is recognized as the 



P 



f 

median of the scale bin with measure p s = s _2 AT 1 , such that f 2 j^ a p s ds = 1/2. The edges for odd 
parity are handled similarly with respect to the half-width bins. For either parity, the bin whose lower 
boundary is a = is the one that causes trouble, as A s = oo in that case. One can contrive to neglect 
that bin for odd parity by subtracting the mean of the signal y(t) — > y(t) — (y(t)) t , but that procedure 
does not repair the problem for even parity. 
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Figure 3. Comparison of the mapping over odd parity axes of frequency / and scale s = f" 1 
of the discrete spectral power of a signal with two sinusoidal components for / G [0, f c ] in 
(a) and (b) and for / G [f c , 2f c ] in (c) and (d) with the central frequencies indicated by x . 
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Table 1. Numerical results from the evaluation of Figure 3. 
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To illustrate the difficulty, in Figure 3 we display the mapping of the spectral power for a signal with 
duration T = 20, a sum of two sinusoids at frequencies 1/40 and 1/4 normalized to unit energy E y = 1. 
In panels (a) and (b) we show the mapping for the one-sided transform of order N = 20 and odd parity 
over domain / G [0, f c ]. While the higher frequency's contribution remains apparent, that of the lower 
frequency has been washed away by the measure factor. The dof carried by the lowest frequency bin are 
unrecoverable from the mapping over scale on account of the infinite bin width. In panels (c) and (d) we 
repeat the procedure but for domain / G [/ c , 2/ c ]. This time, all the dof are properly mapped, thus the 
energy and reconstruction theorems are satisfied. In Table 1 we give some numerical results from the 
evaluation shown in the figure. In short, there is nothing to be gained by working on the scale axis rather 
than frequency other than a headache in dealing with the effect of aliasing, as one requires the same set 
of basis functions and minimal order A^ min to satisfy the fundamental theorems of spectral analysis in 
either case. 



3. Gabor and Morlet Transforms 
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Let us turn now to the consideration of how to transform the signal y = y(t) into a time-frequency 
representation Y = Y(f,t) by modulating the exponential oscillations of the Fourier basis with a 
discrete window function, Q(f,t) — > ty(f,t) = $(t)0(/,t). The Gabor transform is defined by the 
use of a Gaussian window with decay parameter a, which we will notate as oc exp~ 7r / 2 (t 2 /cr 2 ) = 
e ' . Its normalization to unit energy $ $/($|$)J /2 depends upon the duration of the 
window, parametrized by its half-width r. In the continuum, one can write & 2 (t)dt oc a and 
f^ T $ 2 (t)dt oc cr erf {ix 1 I 2 t / o); however, in the discrete setting with index t E [— r, r], the window 
energy Y^ t <3> 2 (t)A t is not expressible in closed form. For the one-sided spectrum, $ — > \/2 $ so that 
its energy equals 2. To fully specify which Gabor transform is being used, one must state the values of 
N, P, a, and r. We will return to the question of finding the minimal order of the Gabor transform for 
a signal of duration T in Section 5; for now, we will assume that N is greater than N min of the previous 
section so that the fundamental theorems are satisfied. 

There are two alternatives for the definition of the phase convention in the transform, which differ in 
whether the phase is expressed relative to the origin of the signal's time axis or that of the window: 

T 

Y(n,t) = H-t')e^[-f n (i + t')}y(i+t')A t , (7) 

t' = -T 
N' T 

§(t) = ^exp^(/„t) Ht')Y(n,t + t')A t A f , (8) 

71=1 t' = -T 

where the time axis for the transform coefficients carries an index t E [1 — r, T + r], or else 

T 

Y(n,t) = H-t')exp^(-f n t')y(i + t')A t , (9) 

t' = -T 
N' T 

kt) = Yl QWexp^i-MYfat + ^AtAf. (10) 

n=l t'=-T 

The latter is more familiar, but either is equally valid in terms of satisfying the energy and reconstruction 
theorems. The sign of the argument to the window function is chosen deliberately so that these 
expressions remain valid for the case of a window which is not symmetric in t' . The signal is zero 
padded for values of t + t' outside the original domain of [1,T], as no other value could be assigned 
without changing the energy hence information content of the signal. For either phase convention, the 
temporal spectral power is defined as P(n } t) = \Y(n, t) \ with units of energy per time per frequency. 

1/2 A 1/2 

(The normalization of the window introduces a factor of A t so that Y carries units of u y u t .) For 
comparison to the signal energy E y , let us evaluate the spectrum's energy as Ey = J2 n Sf? P( n i t)A t A r 

2 

and the reconstruction's energy as Es = \y(t)\ A t . 

For this section and the next, let us consider a particular real signal y{t) comprised of 4 sinusoids with 
unit amplitude and non- stationary frequency of duration T = 200. The variation in the instantaneous 
frequencies is assigned an amplitude of 5% and a period of T. Phases for the oscillations and variations 
are selected randomly. In Figure 4 we show the results of the one-sided Gabor transform with N = 200 
and P = using window parameters of a = 2^pn and r = 12. Panel (a) displays the contours of the 
energy density P(n, i) over axes of time and frequency. Also indicated are the instantaneous frequencies 
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Figure 4. One-sided Gabor transform of a real signal as described in the text. The power 
spectrum is shown in (a) with the signal's instantaneous frequencies indicated by the dashed 
lines. The marginal spectral power is shown in (b) as x , as is the marginal temporal power 
in (c), and each is compared to the convolution of the window and signal energies indicated 
by +. The absolute value of the reconstruction residual is displayed in (d). 




used to generate the data. Panels (b) and (c) display the marginal energy densities over axes of frequency 
and time, respectively. The marginal energy densities are evaluated by taking each sum over an axis 
independently, i.e. P(n) = ^2fP{n, t)A t and P(t) = J2 n P( n i ^)Af- P ane l (d) gives the absolute value 
of the residual r(t) = y(t) — y(t), which is on the order of the machine precision ~ 1CT 16 . The ratios 
Ey j Ey and E* / E Y are equal to unity to the same precision, thus the fundamental theorems are satisfied. 
At these parameter values, the spectral resolution is poor while the temporal resolution is sharp. 

The marginal densities of the Gabor power spectrum can be compared to the convolution of the 
window's energy density with that of the signal in either the temporal or spectral representations. 
Looking first at the temporal representation, one finds the relation 

T 

P(t) = ® 2 (-t')y 2 (t + t')A t , (11) 

t'=-T 

where $ in this case is normalized to unit energy. In the spectral representation, the convolution is taken 
in the modular (periodic) sense, most easily performed using the two-sided transform. Let y(f n ) be the 
two-sided Fourier transform of the signal with order N and parity P, and let &(f' rn ) be the transform 
of the window function, again normalized to unit energy, with the same order and odd parity such that 
M = 2N + 1. Then P{n) in the two-sided sense n E [1, 2N + P] can be written 

M 

P{n) w |<&(M + 1 -m)| 2 |£[mod(n + m - N - 2,2N) + l}\ 2 A f , (12) 

m=l 
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Figure 5. One-sided Morlet transform of a real signal as described in the text. The power 
spectrum is shown in (a) with the signal's instantaneous frequencies indicated by the dashed 
lines. The marginal spectral power is shown in (b) as x , as is the marginal temporal power 
in (c). The absolute value of the reconstruction residual is displayed in (d). 




where mod(a, b) = a mod b and with respect to the half-width of the edge bins, from which the values 
corresponding to the one-sided transform can be extracted and doubled. The relation above is written as 
an approximation because there is some subtlety to the evaluation of the RHS. 

So far, we have imposed no condition on the window duration T = 2r + 1 other than being odd 
for integer r, and indeed, there is none. In the lower limit r = 0, the basis functions are constant, 
\l/(0, n) = y/2 (or 1 for the two-sided transform), and all frequency resolution is lost, yet the fundamental 
theorems are still satisfied for any order N > 1. In the opposite extreme, for example r = T such that 
T" ^> T with a — oo, the basis functions become essentially those of the Fourier transform with minimal 
order iV min . Simply put, the duration of the window in the short-time Fourier transform need not be short. 
When the temporal bandwidth of the window approaches or exceeds that of the signal, the evaluation 
of the RHS in Equation (12) becomes suspect if the order N is not sufficient to resolve both the signal 
and the window; the practical solution is to increase N. Furthermore, one can verify that the Gabor 
transform remains well behaved (if not particularly useful) for the case a — > ia so that the window 
grows exponentially rather than decaying. The only requirement is that the window be real valued and 
normalized explicitly, i.e. discretely, to unit energy, times 2 for the one-sided transform. 

We can now introduce the Morlet transform by promoting the window parameters from constants to 
functions of the central frequency, — >■ $ n (t). Let a n = a/ f n , and with respect to the discrete 
sampling in time, let r n = \t/ f n ~\ for constants a and r. The parameter r itself need not be an integer, 
as r n is what defines the window duration for bin n. A difficulty with r„ arises immediately for odd 
parity when f n = 0; one can ameliorate the situation by using r n from the lowest positive bin in that 
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case. Again, either the phase convention of Equations (7-8) or (9-10) can be used. In Figure 5 we show 
the results of the one-sided Morlet transform of order 200 and even parity with parameters a = y/ir and 
r = 6 chosen so that the window of the previous Gabor transform corresponds to the Morlet window at 
the critical frequency f c . The power spectrum in panel (a) has the marginal densities shown in panels 
(b) and (c); however, this time the ratios E Y /E y and E~/E Y differ from unity on the order of 1%, 
interestingly by nearly the same value. Likewise, the absolute value of the residual displayed in panel 
(d) is on the order of 1%, comparable to the Farge method [2] as reported by Torrence and Compo [1]. 
While the lowest frequency component has been well resolved, the spectral resolution of the upper 
frequencies remains poor. Furthermore, without a single window in operation, one cannot perform the 
comparison of the marginal densities with the convolution of the window and signal energy densities 
corresponding to Equations (11) and (12). 

Let us insert here some comments on what is called the admissibility condition. The statement often 
is made by authors that the wavelet must have a mean of zero so that its Fourier coefficient at / = 
vanishes. While that remark might apply in the continuum, it is not appropriate for discretely sampled 
wavelets of finite duration. If one examines in detail the proof of the admissibility condition [21], one 
finds that it relies crucially on the property of continuity. In full, the admissibility condition states that 
if the wavelet is continuous in time with continuous Fourier transform, then its mean must be zero. 
In the context considered here, neither of those conditions is met: in the temporal representation the 
wavelet is a piece-wise constant function with jump discontinuities at the edges of the temporal bins, 
and likewise for the frequency representation which, while of arbitrary resolution beyond the minimal 
order, must nonetheless be evaluated discretely for any practical analysis of data. If one computes 
the leakage functions for the discrete Morlet basis [22], one finds an artificial discontinuity is induced 
by the imposition of the zero mean condition. As a final argument against subtracting the mean of the 
discrete basis functions, note that the Gabor transform works perfectly well for any values of the window 
parameters without including such procedure. 

There are many suggestions found in the literature for improving the performance of the discretely 
implemented continuous wavelet transform. Let us examine a few of them here, with their absolute 
residuals displayed in Figure 6 and a numerical summary in Table 2. For panel (a) the mean of the 
signal is subtracted before entering the spectral analysis. As one can see, not much changes, but for 
consistency of comparison the remaining panels will also use the mean subtracted signal. Noticing that 
the ratios E Y / E y and E~/E Y are very nearly equal, one can renormalize the spectrum and reconstruction 
a posteriori by their geometric mean C = (E~/ Ey) 1 / 2 , yielding the residual shown in panel (b). This 
renormalization prescription works best when there is not much energy in the upper portion of the 
frequency axis [22]. Perhaps more familiar is the prescription to shift the central frequency 2ir — > u\. 
At these window parameter values, a shift of cui/2rr = 1.0127 is found to work well [23], giving the 
residual shown in panel (c). If one then applies the renormalization prescription, one gets the residual 
displayed in panel (d). While there has been a noticeable improvement of three orders of magnitude, the 
root- mean- square residual remains far above the machine precision. 



4. Layered Window Transform 
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Figure 6. Absolute residuals for several methods of improving the transform response as 
described in the text. In panel (a), the mean of the signal is subtracted; in panel (b), the mean 
is subtracted with renormalization of the coefficients; in panel (c), the mean is subtracted 
with a shift of central frequency; and in panel (d), the mean is subtracted with frequency 
shift and renormalization. 
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Table 2. Numerical results from the evaluation of Figure 6. 



panel 


(a) 


(b) 


(c) 


(d) 




1.8215e-02 


2.6142e-03 


1.8213e-04 


4.1799e-05 




1.2885e-02 


-1.7468e-06 


1.2674e-04 


-4.4658e-10 


E^/E Y - 1 


1.2889e-02 


1.7468e-06 


1.2674e-04 


4.4658e-10 



The primary feature of the Morlet transform is that it uses a different window function $ n for each 
bin along the frequency axis. That is also its main difficulty, as each bin corresponds to a particular 
instance of the Gabor transform, any of which would work fine independently, but when the separate 
bins are isolated and combined, there is no reason to suppose that their marginal sums will equal the 
signal energy. The correspondence between the Morlet and Gabor transforms becomes more apparent 
if one redefines the window decay o n relative to the window width r n rather than the central frequency 
f n , in which case regions along the frequency axis of the Morlet spectrum are drawn from a particular 
Gabor spectrum. Simply pasting together pieces from different Gabor transforms from the outset appears 
doomed to failure, and it is surprising that the Morlet transform in the discrete setting works as well as it 
does. While we have not given up entirely on the Morlet transform, if one wants to make quantitative use 
of the results of a spectral analysis, the energy and reconstruction theorems must be satisfied precisely. 
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How, then, are we to devise a multiresolution analysis that satisfies the fundamental theorems of 
spectral analysis? The answer lies in focusing our attention on the distribution of energy in the window 
function. Suppose we have one unit of energy, an "atom" so to speak, distributed with temporal power 
&\{t) in the continuum, which we wish to use as the basis for a multiresolution analysis. Let the window 
at scale r be defined by $ T (£) oc $yi(t/r) for positive integer r, which scales the relative unit of time by 
integral amounts. One can supplement this definition for r = by letting $o(0) = 1 and otherwise. 
For an exponential decay with scale invariant parameter a, one can write & T (t) oc exp~ 7r / 2 (t 2 /V 2 cr 2 ), 
where the normalization to unit energy is done explicitly over the integers t E [— r, t\. (Alternately, one 
could hold the window width fixed for all r.) The layered window $l(^) is then defined to be the sum 
in terms of energy, not amplitude, of the windows & T (t) over some set of scales r £ [n, tl] with L 
members, normalized to unit energy, $|(t) = L^ 1 J2 T F° r a one-sided transform, $ L — > \[2^ L 

of course. This single window, with duration T' = 2t l + 1, is to be used for all the frequency bins in the 
windowed Fourier transform. 

Since the atomic window $^4 is arbitrary, the distinction between summing the window energies $ 2 
versus summing the window amplitudes $ T is really just a matter of interpretation. However one defines 
the sum over scales, the layered window eventually is normalized explicitly to one unit of energy, or two 
units for a one-sided transform. Using the construction above in terms of energy, the amplitude of the 
layered window is the root-mean- square of the amplitudes of the scaled windows, while the construction 
summing the window amplitudes is not so simply normalized. With respect to its physical motivation, 
energy is the quantity that ultimately gets quantized, not amplitude, which is a property reflected in our 
favored definition of the layered window. The selection by the investigator of the set of scales {r} to use 
in the construction specifies the domain of energy distributions to which the transform will be sensitive. 

Let us return now to consideration of our four frequency test signal, with its mean restored. For 
comparison with the previous Gabor and Morlet transforms, let the decay parameter here be a = y/n/G, 
and let r G [12, 96]. The results of such analysis are displayed in Figure 7. Compared to the Morlet 
transform, the spectral resolution is much improved in the upper portion of the frequency axis, as seen 
in panel (a). Since only one window is in operation, the marginal densities of panels (b) and (c) can be 
compared to the convolution of the window and signal energy densities in either the temporal or spectral 
representations. The absolute residual shown in panel (d) is on the order of the machine precision, with 
maybe a slight accumulation of truncation errors. The ratios Ey/E y and E^jEy equal unity to machine 
precision. The layered window transform, by virtue of using a single window for all frequencies, satisfies 
the energy and reconstruction theorems to the accuracy of one's computational device. 

As a final comparison, let us put the Fourier power spectrum of the signal alongside the marginal 
spectral energy densities of the Gabor, Morlet, and layered window transforms, displayed in Figure 8. 
The Fourier spectrum in panel (a) is the typical mess one gets when analyzing non- stationary signals, 
but its net energy does equal the signal energy. The Gabor spectrum in panel (b) has lost its frequency 
resolution on account of the tight temporal bandwidth of the selected window parameters yet retains 
the correct normalization. The Morlet spectrum in panel (c) resolves the low frequency portion of 
the spectrum but not the high frequency region and does not have the correct normalization. The 
layered window spectrum in panel (d) is normalized to the signal energy and resolves both the low 
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Figure 7. One-sided layered window transform of a real signal as described in the text. The 
power spectrum is shown in (a) with the signal's instantaneous frequencies indicated by the 
dashed lines. The marginal spectral power is shown in (b) as x, as is the marginal temporal 
power in (c), and each is compared to the convolution of the window and signal energies 
indicated by +. The absolute value of the reconstruction residual is displayed in (d). 




and high frequency parts of the spectrum while maintaining sufficient temporal resolution to be useful 
in identifying non- stationary features in the signal. 



5. Irregular Sampling and Minimal Order 

Let us now consider the case of an irregularly sampled signal y(t), with regard to the analyses 
by Lomb [24] and Scargle [25]. Irregular sampling has long vexed practitioners of wavelet analysis, 
leading to a variety of suggestions for how to overcome its difficulties. Because of the nature of stellar 
observations, the astrophysical community often presents data which contains gaps or is otherwise on an 
irregular time axis. The method by Foster [26] focuses on the normalization of the gapped wavelet basis 
functions, while the method by Frick et al. [27] focuses on the admissibility condition; taken together, 
those ideas lead to the edge adapted algorithm proposed by Johnson [28]. More familiar perhaps is 
the lifting scheme by Sweldens [29], which operates in the context of the dyadic wavelet transform 
implemented in terms of finite impulse response digital filter banks. 

For simplicity, we will assume that the samples each have a uniform duration A t which is not greater 
than any of the intermeasurement periods, so that the time metric D t remains proportional to the identity 
matrix; otherwise, a suitably generalized D t must be used. Let the observation times be given by the 
vector t = t d indexed by d G [1,-D], where D is the total number of samples, in some unit u t not 
necessarily equal to A t such that the values td are integers; u t is the resolution of the time measuring 
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Figure 8. Comparison of the Fourier power spectrum in (a) to the marginal spectral energy 
densities of the Gabor transform in (b), the Morlet transform in (c), and the layered window 
transform in (d). 
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apparatus, thus A t must also be an integer. Similarly, the measurements can be written as the vector 
y = y d so that the signal energy remains E y = y^D t y in units of u^u t . Missing values, i.e. measurements 
at integer times not in t, are effectively treated as zero. Let us look first at how the discretized Fourier 
transform is modified in this case. 

To identify the Nyquist critical frequency, one must find the lowest positive Fourier basis function 
which is entirely real (technically up to a constant phase) over the given set of observation times [30]. 
That procedure is easily accomplished by shifting and scaling the time axis t — > t' such that f c = 
(2u t /) _1 in the new time units. If one defines t g to be the greatest common divisor over the set of 
intermeasurement periods, the new time axis can be written t' d = (t d — t\)/t g + 1 in units of Uf = t g u t 
such that t' d G [1,T] remains integer valued. At the critical frequency, the Fourier basis function is 
exp m (2f c t' d ) = ±1, and the Fourier spectrum has period 1 in units Uf = u^, 1 . Having found the units 
in which f c = 1/2, let us drop the prime distinguishing the scaled time axis in the following, and for 
convenience let us suppose A t — 1 in the scaled units; otherwise one must be extra careful with the 
scaling of the energy units. 

The next task is to determine the minimal order for satisfaction of the energy and reconstruction 
theorems. Let Nr, = \D/2], and let Nt = \T/2]. For N < Nn there are an insufficient number of dof 
in the discrete Fourier transform to fully represent the information content of the (real) signal, thus the 
fundamental theorems are not satisfied in general. Likewise, for iV > N T there are more than enough 
dof to represent the signal; Nt gives the minimal order of the analogous regularly sampled signal where 
the missing entries are assigned the value zero. For N D < N < N T , the Fourier transform might work, 
depending upon the particulars of the irregular sampling. To satisfy the reconstruction theorem y = y 
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Table 3. Comparison of the norm of the residual of the Fourier transform to the distance 
from the quality matrix Q to the identity \ D at various orders N and even parity for an 
irregularly sampled signal with D = 10 and T = 28 such that Np = 5 and Nt = 14. 



N 


5 


6 


7 


8 


9 


HQ 


\t\\f 
— \d\\f 


3.5943e+00 
2.4495e+00 


5.1760e+00 
2.8284e+00 


3.0368e+00 
2.4495e+00 


5.5110e-15 
5.6505e-15 


3.0148e+00 
2.0000e+00 


N 


10 


11 


12 


13 


14 


HQ 


M\f 
— \d\\f 


6.6842e-01 
1.4142e+00 


5.4143e-15 
4.3648e-15 


2.6771e+00 
2.0000e+00 


6.241 5e- 15 
4.2983e-15 


7.4742e-15 
4.9648e-15 



for y = Qy in the one-sided case, one must have Q = Re0D/0^D t = 1^ to the order of machine 
precision; the energy theorem also is satisfied when that condition is met. (For a complex signal, the 
entire matrix, not just its real part, must be utilized, and of course D f must cover one full period of the 
frequency axis.) The basis functions are evaluated only at the observation times Q(d, n) = exp l27T (f n td), 
because those are the only locations at which the reconstruction theorem must be satisfied. In Table 3 we 
compare the norm of the reconstruction residual r = y — y for some signal with D — 10 and T = 28 to 
the distance between the quality matrix Q and the identity I D according to the Frobenius metric at various 
orders N with P = 0. To produce a smooth picture of the spectral content, one should take N 3> N T . 

Turning now to consideration of the Gabor transform (and by extension all windowed Fourier 
transforms with a fixed $), the construction of the quality matrix is a bit more complicated, owing to the 
convolutions along the time axis. It is commonly remarked that the number of dof along the frequency 
axis is multiplied by the number of times the window duration T" = 2r + 1 fits within the signal duration 
T, but the situation is more subtle than that, as one must also account for the bandwidth of the window 
given by the decay a. Beginning with the case of regular sampling D = T, using the phase convention 
of Equations (9-10), one can write the composition of the inverse and forward transforms as 

N' T T 

§(t) = Re J2 J2 exp^H/Wt') E exp i27r (-fJ")y(t + t' + t")A t AAj(tta) 

n=l t'=-r t"=-r 
2t N' 

= Re E E wm-t")Y. e ^ 2 ^-frf')y(t+n^f^At (i3b) 

t»'=-2r t'+t"=t"' n=l 
2t 

= q(t"')y(t + t"')A t , (i3c) 

i"' = -2r 

which yields one row in the quality matrix Q indexed by t'" E [— 2r, 2r] indicating the diagonal where 
q(f") appears such that Q has the form of a Toeplitz matrix of rank T. For a symmetric window $(—£') = 
$(t') written as a diagonal matrix <t» such that W = <t>0, the values q(t"') can be extracted from U^D/U't . 
For an irregularly sampled signal T > D, one retains only those rows and columns of Q corresponding 
to actual data entries, Q — > Q(t, t). In Table 4 we compare the residual norm of the Gabor transform 
with r = 12 and two values of a to the distance from Q to Id at various N with P = for an irregularly 
sampled signal of duration T = 100 with D = 30 entries. As N approaches N D from below, the quality 
of reconstruction improves until the residual is on the order of machine precision. 
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Table 4. Comparison of the norm of the residual of the Gabor transform with width r = 12 
and decay a indicated to the distance from the quality matrix Q to the identity I D at various 
orders N and even parity for an irregularly sampled signal with D = 30 and T = 100 such 
that N D = 15 and N T = 50. 
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|r|| F 
-Id 


If 


5.4222e-04 

J.iJjOC Ul 


2.8233e-05 


5.1908e-07 
4 4Qfi?e-ft7 


1.1113e-08 


7.3918e-ll 

5 1fi10e-1 1 

J.1U1 1 1 


2^ 


N 


11 


12 


13 


14 


15 




||Q 


Hf 
-Id 


F 


2.7305e-13 
1.6628e-13 


1.2477e-14 
4.3047e-16 


6.2789e-15 
4.1495e-16 


5.8445e-15 
2.8105e-16 


7.2210e-15 
2.2323e-16 


a 


N 


6 


7 


8 


9 


10 




HQ 


Hf 
-Id 


F 


4.5343e-01 
4.3784e-01 


2.6201e-01 
1.8839e-01 


7.5299e-02 
6.5223e-02 


3.4864e-02 
2.2263e-02 


6.9207e-03 
4.8319e-03 




N 


11 


12 


13 


14 


15 




HQ 


Ml* 
-Id 


F 


1.0987e-03 
6.6716e-04 


7.8620e-05 
6.9627e-05 


1.0067e-14 
6.0970e-16 


1.2076e-14 
7.7092e-16 


1.2964e-14 
1.0855e-15 



Table 5. Comparison of the norm of the residual of the Morlet transform with width r = 6 
and decay a = y/n to the distance from the quality matrix Q to the identity I D at various 
orders N and even parity for a regularly sampled signal with D = T = 30 such that N D = 

N T = 15. 



N 


6 


7 


8 


9 


10 


iiq 


— 'd||f 


3.347 le+00 
2.9527e+00 


2.9397e+00 
2.4686e+00 


2.7428e+00 
2.0993e+00 


2.2582e+00 
1.8244e+00 


1.9475e+00 
1.5788e+00 


N 


11 


12 


13 


14 


15 


HQ 


— 'd||f 


1.791 le+00 
1.345 le+00 


1.5917e+00 
1.1139e+00 


1.3310e+00 
8.7648e-01 


9.3137e-01 
6.3175e-01 


5.3677e-01 
4.0199e-01 



For the Morlet basis functions, the construction of v|/ is more involved, as ^f(t', n) = $(£', n)Q(t', n). 
Nonetheless, one can evaluate Q from U^D/U't for any rank T. As a function of order N, one finds 
that the deviation ||Q — \t\\f decreases until N > T/2, after which it bottoms out on the order of a few 
percent times T. In Table 5 we compare the norm of the residual for the Morlet basis using a = ^fix 
and r = 6 for some regularly sampled signal with D = T = 30 to the deviation of Q from 1^ at 
various N with P = 0. If one wishes to investigate the feasibility of devising a Morlet basis with perfect 
reconstruction, the evaluation of Q is where to start. 

Let us close this section by looking at the analysis of the signal from the previous section but with 
the number of measurements reduced by a factor of a third (N D = 134) for the same duration T = 200. 
The evaluation of the layered window Fourier transform proceeds as before, with the understanding that 
values of y(t + t') at times not in t are treated as zero; the time axis for the spectrum Y(n, i) includes 
every integer i E [1 — t l , T + r L ]. Using the same parameters as before, a = ^/tx '/6 and r G [12, 96] 



Version May 23, 2013 submitted to Axioms 



17 of 24 



Figure 9. One-sided layered window transform of an irregularly sampled signal as described 
in the text. The power spectrum is shown in (a) with the signal's instantaneous frequencies 
indicated by the dashed lines. The marginal spectral power is shown in (b) as x, as is the 
marginal temporal power in (c), and each is compared to the convolution of the window and 
signal energies indicated by +. The absolute value of the reconstruction residual is displayed 
in (d). 




with N = 200 and P = 0, in Figure 9 we display the power spectrum, its marginal densities, and the 
reconstruction residual for our irregularly sampled signal. While the loss of information has affected 
the appearance of the spectrum relative to the regularly sampled case, it nonetheless remains a faithful 
representation of the information which is available. By virtue of using a single fixed window & L , the 
layered window Fourier transform satisfies the fundamental theorems of spectral analysis even when 
there are gaps in the observation record. 

6. Window Comparison 

Let us conclude the analysis with a comparison of the temporal and spectral bandwidths for several 
types of window function, as well as their estimates of the power spectral density carried by a test signal 
with a little more complexity than we had before. The test signal is regularly sampled with T = 200 and 
four component frequencies, but now only two components are present in the first half of the duration, 
while the other two are in the second half. An abrupt transition occurs at the midpoint of the duration. 
The amplitude of the frequency variation is now 10%, with a period of T/2, so that each component 
covers a broad range of instantaneous frequency. The order for all the transforms in this section will be 
the minimal Fourier order of the signal N = 100 with even parity P = 0. 
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Figure 10. Comparison of the temporal and spectral energy densities of the various windows 
arranged according to Table 6. 
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Table 6. Comparison of the temporal and spectral bandwidths for various windows, as well 
as their rms reconstruction residual for the test signal as described in the text. 



Window 


a 


T 


h 


t 2 


fi 


/ 2 


4^2/2 


\\*\\f 


Gabor 


2^ 


12 


-0.000 


1.414 


-0.000 


0.056 


1.000 


0.000 


Gabor 




96 


-0.000 


11.314 


0.000 


0.007 


1.000 


0.000 


Layered 


A/6 


[12,96] 


0.000 


6.990 


0.000 


0.014 


1.244 


0.000 


Random 




24 


-2.836 


13.718 


0.000 


0.130 


22.442 


0.000 



The parameters for the four types of window considered are displayed in Table 6, comprising of two 
Gabor windows at either end of the range in scale for the layered window also considered, in addition to 
a random window with an arbitrary duration. The first moments of time and frequency are indicated by t\ 
and fi respectively, and the bandwidths (square root of the second moment about the first moment) by t% 
and f 2 . The two Gabor windows are observed to minimize the Fourier uncertainty relation £2/2 > 1/47T, 
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Figure 11. Comparison of the power spectral density of the test signal described in the text 
using the various windows arranged according to Table 6. The instantaneous frequencies 
used to generate the signal are indicated by the dashed lines. 




t t 



thus in that sense are optimal, while the layered window has a slightly larger bandwidth product. The 
random window has a bandwidth product which is considerably larger, yet all these windows produce 
a valid discrete transform pair that satisfies the fundamental theorems of spectral analysis as indicated 
by the residual of reconstruction ||r|| F . The energy densities in the temporal and spectral representation 
used to evaluate the bandwidth products are shown in Figure 10, where one can see the trade-off in 
time/frequency resolution in action. 

The power spectral density estimates evaluated from the test signal of this section using the various 
windows are displayed in Figure 11, arranged from (a) to (d) according to Table 6. The power spectrum 
given by the layered window in (c) does indeed incorporate features of either Gabor transform shown 
in (a) and (b); in a sense, it has combined the Gabor transforms over its range of scale in a way that 
preserves the energy and reconstruction theorems. With some algebraic manipulation of the expressions 
defining the layered window Fourier transform one might be able to show that explicitly, but for now that 
statement is intuitive speculation. Interestingly, the spectrum given by the random window is not much 
different from the others, driving home the point that the shape of the window truly is arbitrary as long 
as it is normalized appropriately. 

7. Discussion 

The primary result of this investigation is that the windowed Fourier transform has far more flexibility 
and utility than it is usually credited with. To achieve satisfaction of the energy and reconstruction 
theorems in the discrete setting, the only requirements on the window are that it be real, nonnegative, 
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and normalized explicitly to unit energy, or two units for a one-sided transform. While we have looked 
only at symmetric windows here, one can easily verify that the expressions for the transform pair, either 
Equations (7-8) or (9-10), remain valid for a window that is not symmetric in t'; with attention to the 
details of relocating the temporal bins, the coefficients can be assigned to the time corresponding to 
the peak of the window rather than its midpoint. Similarly, a window with an even length duration T' 
could be used if one is careful with the definition of the phase and the location of the bins. In fact, 
any nonnegative function can be used for <E>, or even a list of arbitrary numbers, as long as it is suitably 
normalized. 

The flexibility in $ allows one to define a multiresolution spectral analysis in terms of an atomic unit 
of energy $^(t/r) evaluated over a range of scales r. These multiple scalings of &a are combined into 
a single layered window $l which is applied to the entire frequency axis; the phase component of the 
basis is independent of the window function. Because only a single window is used, the fundamental 
theorems of spectral analysis are satisfied. With a Gaussian the form of $l very closely resembles 
a Lorentzian function, which expresses uniformity over scale when expressed as an angle. The shape of 
$ l can be tailored to the needs of the investigation through inspection of its leakage function. Similarly 
to the wavelet transform, the trade-off between resolution in time or frequency is under the control of the 
investigator. Any structure seen in the spectral coefficients of the signal is understood to be conditioned 
on the selection of the window function; the answer one gets depends upon how the question is asked. 

The difficulties faced by the continuous wavelet transform exemplified by the Morlet basis, with 
respect to the fundamental theorems of spectral analysis, can be summarized in the evaluation of its 
quality matrix Q. Perfect reconstruction in the discrete setting requires Q to equal the identity matrix to 
the precision of one's computational device. The various adjustments looked at above which improve 
the reconstruction residual must be bringing Q closer to that form; however, none of them achieve the 
required precision. For an integral transform pair to have quantitative significance, one must demonstrate 
that the distance from its Q to I is numerically zero. This investigation has proposed an alternate approach 
to multiresolution analysis which does satisfy the energy and reconstruction theorems while improving 
upon the resolution properties of the Gabor transform. 

The effects of aliasing and of irregular sampling are easily understood with regard to the periodicity 
induced in the Fourier spectrum. The Nyquist interval is given simply by the span between frequencies 
which are indistinguishable over the stated measurement times; assigning a value of unity to that range 
such that f c = 1/2 is just a matter of choosing the appropriate units for time. The implementation of 
the windowed Fourier transform is unaffected beyond zero padding the signal at locations of missing 
measurements, whose justification is that no other procedure leaves the energy, hence information 
content, unchanged. One should observe that the ends of a regularly sampled signal are treated the 
same way. On that note, we can remark that there is no cone of influence for the windowed Fourier 
transform, as every spectral coefficient is important to the satisfaction of the fundamental theorems, even 
those outside the temporal domain of the data. 

The minimal order iV min required for a faithful representation of the signal depends upon the temporal 
bandwidth of the window and the duration of the data. Pinning that number down for the windowed 
Fourier transform is a bit tricky, and if needed is best evaluated explicitly through inspection of Q as 
a function of order N. What one can say in general is that 1 < iV min < N T , where the bounds are 
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Figure 12. Comparison of the lowest ^f(t', 1) and highest ^(f, N) frequency basis functions 
at order N = 50 with P = using the parameters found in the text, with the real part 
indicated by x and the imaginary part by +. The Fourier basis is in (a) and (b), the Gabor 
basis is in (c) and (d), the Morlet basis is in (e) and (f), and the layered window basis is in 
(g)and (h). 
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given by the minimal orders for the fully temporal and fully spectral representations of the regularly 
sampled signal, respectively. Irregular sampling complicates matters by introducing an order N D < N T 
such that Nd < N min < Nt for the full-length Fourier transform and A^ min < Np for the windowed 
Fourier transform. For most practical cases of data analysis, one is interested in producing a smooth 
plot of the spectral content, thus one would use an order N N min . To verify the marginal density 
of the windowed power spectrum in comparison to the convolution of the window and signal energy 
densities in the spectral representation, one requires an order N sufficient to resolve both the signal and 
the window temporal durations. 

Thus far we have heard nary a peep from the actual basis functions used in this analysis, so let us 
close the discussion by looking at some of the stars of the show. In Figure 12 we display the real and 
imaginary parts of the basis functions, normalized to unit energy, for each of the transforms considered 
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above at order N = 50 with even parity P = at the lowest and highest entries on the frequency axis, 
i.e. fi = 1/200 and /at = 99/200. The Fourier transform has no intrinsic window, so it is shown for 
a duration T 1 = 101 which spans T < T'. The window parameters chosen are those from the previous 
sections: o = 2y/n and r = 12 for the Gabor transform, a = y/n and r = 6 for the Morlet transform, 
and a = y/n/6 and r 6 [12, 96] for the layered window transform. One can observe that the transforms 
whose Q equals the identity share the property that these basis functions are phase analogues over the 
discrete sample times, whereas the Morlet basis functions are scale analogues. By phase analogue we 
mean if 1) = a x + ib x , then ^(t', N) = ai(-l)*' + ib 1 (-l) t ' +l , a feature shared by all positive 
frequency pairs whose midpoint is / = 1/4; the Morlet basis differs from the others in this regard. 
Whether that property is required for satisfaction of the fundamental theorems is unclear, but the results 
of this analysis suggest that it is. 

8. Conclusions 

In this article we have compared the Fourier, Gabor, and Morlet transforms in the fully discrete 
setting applicable to the analysis of sampled data. While the Fourier and Gabor transforms satisfy the 
fundamental theorems of energy conservation and perfect reconstruction, the Morlet transform does not. 
The magnitude of the residual can be related to the distance from the quality matrix to the identity in 
all cases, such that the minimal order for reconstruction can be determined for the Fourier and Gabor 
transforms. Various methods of improving the response of the Morlet transform are considered; however, 
none of them achieve the desired precision on par with the truncation error of one's computational device. 

An alternate approach to multiresolution analysis is proposed which constructs a single layered 
window from multiple scalings of some atomic unit of energy. This layered window transform satisfies 
the fundamental theorems of spectral analysis while providing sufficient temporal resolution to identify 
non-stationary features in the signal. The trade-off between time and frequency resolution is under 
the control of the investigator through the selection of the atomic window and the scales over which 
it is evaluated. The power spectrum of the layered window transform is similar to that of the Morlet 
transform but provides much better frequency resolution. 

The premise behind the wavelet transform is that the low frequency elements of a signal should have 
a much longer duration than the high frequency elements. There are, however, many examples of real 
world signals for which the converse is true. Consider, for example, the digital recording of a kick drum 
and cymbal rhythm such that the low frequency bursts have a relatively short duration compared to the 
ringing at high frequency. For that type of signal, the wavelet transform is going to provide a poor choice 
of basis even if it satisfied the fundamental theorems. The flexibility in the choice of window function 
used in the windowed Fourier transform allows one to tailor its response to the needs of the analysis far 
better. The result one gets depends upon how one defines the resolution of the window, any of which 
provide a valid spectral representation of the signal. 

To make quantitative use of the spectral analysis of some discretely sampled signal, one must 
demonstrate that the fundamental theorems are satisfied. The windowed Fourier transform can be shown 
to satisfy those requirements for any real valued window function that is suitably normalized. The 
potential for gaps in the data, or some other form of irregular sampling, is found to pose no problem 
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once the correct Nyquist interval is identified. Consequently, the windowed Fourier transform can be 
applied to data from a wide variety of sources, such as astronomical observations, which are limited 
physically to an irregular set of observation times, as well as the common case of regular sampling. The 
implementation of a multiresolution spectral analysis in the discrete setting appears to be possible only 
when the multiple scalings of the energy distribution are applied evenly across the frequency axis by 
combining them into a single layered window. 
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